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Elastic capsules in shear flow: 
g '■ Analytical solutions for constant and time-dependent shear rates 
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Abstract. We investigate the dynamics of microcapsules in linear shear flow within a reduced model with 
two degrees of freedom. In previous work for steady shear flow, the dynamic phases of this model, i.e. 
swinging, tumbling and intermittent behaviour, have been identified using numerical methods. In this 
paper, we integrate the equations of motion in the quasi-spherical limit analytically for time-constant and 
time-dependent shear flow using matched asymptotic expansions. Using this method, we find analytical 
expressions for the mean tumbling rate in general time-dependent shear flow. The capsule dynamics is 
studied in more detail when the inverse shear rate is harmonically modulated around a constant mean value 
for which a dynamic phase diagram is constructed. By a judicious choice of both modulation frequency and 
phase, tumbling motion can be induced even if the mean shear rate corresponds to the swinging regime. 
We derive expressions for the amplitude and width of the resonance peaks as a function of the modulation 
frequency. 

PACS. 87.16.D- Membranes, bilayers, and vesicles - 47.15.G- Low- Reynolds-number (creeping) flows 



1 Introduction 

The dynamic motion of soft objects such as elastic mi- 
crocapsules in shear flow represents a long-standing prob- 
lem in science and engineering. It has received increas- 
ing interest recently, in particular due to its relevance to 
biological, medicinal and microfluidic applications. This 
problem is challenging from a theoretical point of view, 
because the shape of these objects is not given a pri- 
ori, but determined dynamically from a balance of h> 
terfacial forces with fluid stresses. New insight has been 
gained due to a plethora of experimental [Jj, 0, d, 0, [j| , 
theoretical [j, E, I I , H I , jicj - fill . [H, [H, [Jj|], and numerical 
[U, M, E3, El Hfla El methods. 

Perhaps the most well-known dynamic state of initially 
spherical elastic microcapsules in shear flow is the tank- 
treading motion also present in fluid vesicles [12, HU HJ, 
[H Sfl H3, El HI, IS Ell) as reviewed in the first two 
chapters of [33]. In contrast to fluid vesicles, microcap- 
sules exhibit a finite shear elasticity, since their membrane 
is chemically or physically cross-linked. This holds both 
for artificial polymerised capsules @ and red blood cells 
(RBCs), whose membrane is composed of an incompress- 
ible lipid bilayer underlined by a thin elastic cytoskeleton 
[33j . For a short time, viscous fluid vesicles can also resist 
shear. 



The resistance to shear leads to qualitatively different 
behaviour, such as preventing the prolate to oblate shape 
transition of viscous fluid vesicles [l7| • Perhaps most sur- 
prisingly, it also leads to qualitatively different instabili- 
ties like wrinkling first observed on polymerised capsules 
0, [l3[ and later as a transient on viscous vesicles [HI . 

When the unstressed initial shape of the cell is not 
spherical, material elements of the membrane are deformed 
when displaced from their initial position. This shape mem- 
ory, suggested for RBCs in Ref. [3j, leads to an oscillation 
of the inclination angle superimposed on the tank-treading 
motion, called swinging, and an intermittent regime be- 
tween tank-treading and tumbling [H, [Til ] . The swinging 
motion of RBCs was studied numerically in Ref. [lj| us- 
ing a boundary integral formulation of the hydrodynamics. 
Later, more comprehensive studies of all dynamic phases 
were performed using both a spectral numerical method 
flcl and an immersed boundary lattice Boltzmann method 
[lg . [2fjj . The phase diagram constructed in Ref. [l8:] ba- 
sically confirmed the qualitative correctness of a reduced 
model [IH at low to moderate viscosity ratios. However, 
both Refs. and [2(J independently contested the in- 
termittent regime at large viscosity ratios as an artifact 
of the reduced model. Instead, in these works the tum- 
bling motion was found to be a transient towards a stable 
swinging motion. Numerical studies of elongated capsules 



2 



Steffen Kessler et al.: Elastic capsules in shear flow 



in extensional flow at high flow rates reveal a novel bi- 
furcation between a spindled and a cusped capsule shape 



induced by compressive stresses [21 



New phenomena are expected when the driving shear 
flow is no longer constant in time and space. Indeed, a 
transient wrinkling phenomenon was observed for fluid 
vesicles in suddenly reversed elongational shear flow [34| , 
where the stress becomes momentarily compressive. For 
spatially varying shear flow produced by structured mi- 
crochannels, a transition from prolate shape to bullet-like 
shape as well as symmetry breaking transitions were ob- 
served in vesicles both experimentally and in simulations 

So far, microcapsule dynamics has only been studied 
in steady shear flow. It is the aim of this paper to in- 
vestigate the effects of modulating the shear rate on the 
dynamics of capsules. In particular, we want to focus on 
the question whether dynamic phase transitions can be in- 
duced by small amplitude oscillations around a fixed mean 
shear rate. Since solving the exact equations of motion nu- 
merically is computationally prohibitive, we constrain our 
investigations to the reduced model of Skotheim et al. [T3 |. 
which will allow analytical solutions in the quasi-spherical 
limit. 

This paper is structured as follows: After reviewing the 
reduced model we derive non-dimensional equations 
of motion valid in the quasi-spherical limit in section [2] A 
numerical study of the dynamics for constant shear rate 
yields the dynamic phase diagram. In the quasi-spherical 
limit the equations of motion can be integrated exactly for 
time-constant shear flow, which is done in section [3] Inte- 
gration is even possible for general time-dependent shear 
flow, which is studied in section|3]in detail. Specialising to 
harmonic modulations of the shear rate around a constant 
mean value reveals that tumbling motion can be induced 
from the swinging regime at certain resonance frequencies 
of the modulation. Both the resonance behaviour and the 
smooth off-resonance background are studied in detail. A 
dynamic phase diagram as a function of modulation am- 
plitude and frequency is constructed. The more intricate 
details of the calculations are shown in the Appendices. 



2 Reduced model 

2.1 Equation of motion 

We investigate a reduced model of an elastic capsule with 
fixed ellipsoidal shape (semi-axes oi, 02, 03, and volume 
V) in a linear shear flow with velocity 



v(x) = jx 2 ei . 



(1) 



depending upon position x = J2i x i e i and shear rate 7. 
The two axes a± and 02 lie in the shear plane with a\ > a 2 
(see Fig. [IJ). Thus the axis 03 is perpendicular to the shear 
plane. The long axis a\ is oriented with inclination angle 
with respect to the direction of the shear flow. The in- 
ner and outer flow have viscosities rf and rf , respectively. 
The membrane can tank-tread with respect to the fixed 




unstressed position. 
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Fig. 1. 2d-cut of a 3d-ellipsoidal capsule with semi-axes a; and 
volume V orientated in an external linear shear flow with outer 
viscosity rf encapsulating a fluid with inner viscosity rf . The 
inclination angle measures the angle between the direction 
of the long axis and the direction of the shear flow. The phase 
angle (f> measures the tank-treading motion. 



ellipsoidal shape, measured by the phase angle 6. This is 
the Keller-Skalak model for a Red Blood Cell Abkar- 
ian et al. [H and Skotheim and Secomb [l4| add an elastic 
energy term (Eq sin 2 cf>) which is due to the tank-treading 
motion and the shape memory effect [3]. Abkarian et al. 4] 
also consider a viscosity of the membrane which effec- 
tively changes the inner viscosity. A Keller-Skalak-type 
Q derivation, which consists of a balance of torque and 
energy, yields the equations of motion for the angles 9 and 

M- 



d t 9 

d t <f> 



2 
2 



2axa 2 7«i 
2-d t (j> +77 — 



2a\ 



■ cos 26 , 



7/3 



E 



■ sin 2(f) — cos 26 



(2) 
(3) 



As the underlying equations of motion are overdamped, 
these equations hold also for time-dependent shear rate 
7 — 7(*)- The geometrical quantities fi depend only upon 
the semi-axes as given explicitly in Appendix rAJ 

In the equations of motion and ([3]) there are seven 
independent parameters, namely the hydrodynamic pa- 
rameters 7, rf, 77°, an clastic parameter Eq and the geo- 
metric parameters ax, 0,2, 03 which determine V, fx, f 2 , 
fs. Three of them can be used to introduce independent 
scales. The volume V defines a length scale, the shear rate 
7 of the external flow a time scale, and the elastic energy 
Eq an energy scale. 

Four independent parameters remain, i.e. the viscosity 
contrast rf /r]° the ratio between hydrodynamic and elas- 
tic energy (^"V^/Eq, the ratio of the short to the long 
axis in the shear plane 02/01, and the ratio of the axis per- 
pendicular to the shear plane to the long axis 03/01. The 
first three can be used to define three equivalent nondi- 
mensional parameters, the shifted nondimcnsional viscos- 
ity contrast 



h rf -h 
-2/3 rf -2/3 



(4) 
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8/tt 

Fig. 2. a) Phase space T = S 1 x S 1 of angles (0,</>). b) By 
cutting the torus along the coordinate axes, it can be mapped 
onto the square [0, ir] 2 on the plane by identifying the pair of 
lines 9 = 0, 9 = n and (f> = 0, <f> — n - a ) an d b) Coordinates 
9, (j>, new coordinates S, A denned by eqns. (|18lll9fl . and two 
typical trajectories (solid lines). Arrows in b) denote direction 
in time. 



(note that /i > 0, /2,/3 < 0, and A > for physical 
values), the capillary number 



and the eccentricity parameter 



arctan 



at 



(5) 



(6) 



Here, a — > corresponds to the spherical case (02/01 — > 
1), and a — > 7r/2 corresponds to the case 02/01 — > 0. It is 
convenient to introduce a dimensionless time r by 



^2 

A 



df . 



(7) 



which can also be done in the case of a positive time- 
dependent shear rate 7 = j(t) > 0. We finally arrive at the 
nondimcnsional reformulation of the equations of motion 



d T 9 
d T cj> 



- cos a d T <p - 
■{x' 1 sin 2^ 



A(l- sin a cos 29), (8) 
f cos 26»), (9) 




Fig. 3. a) Tumbling: The capsule flips with montonously 
changing inclination angle 9 and an oscillating phase angle <j>. 
b) Swinging: The inclination angle 9 oscillates around a con- 
stant value while the phase angle is changing monotonously, c) 
Intermittent motion: Both angles grow without bounds leading 
to a mixture of tumbling and swinging motion, d) Trajectories 
for swinging (black), tumbling (dark grey), and intermittent 
motion (light grey). 



where = x _1 ( T ) > can be time-dependent. The 
phase space is the torus (8, 4>) e T = S 1 x S 1 (each angle 
with period tt, see Fig. [2|. Note that the fourth parame- 
ter 03/01 does not enter the nondimensional equations of 
motion explicitly. 



2.2 Mean tumbling rate and phase diagram 

The solutions of the equations (|8|) and j9|) in the case of a 
time-constant shear flow have been examined in Ref. [l4[ ■ 
The phase diagram consists of three different regimes de- 
pending upon the value of three parameters A, x" 1 , and 
a. As shown in Fig. [3l the capsule can either tumble with 
a monotonously changing inclination angle 9 and an oscil- 
lating phase angle <f> or tank-tread with an oscillating incli- 
nation angle 9 and a monotonously changing phase angle 
<p, a motion called swinging. There is a third regime in 
between, where the capsule both tumbles and tank-treads 
(either successivly or simultanously) , which is called "in- 
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Fig. 4. Typical phase diagrams for different eccentricity parameters a — O.OIO^, 0.333^-, 0.667^ with the dimensionless inverse 
capillary number x _1 an d the dimensionless viscosity ratio A as axes. The mode of motion for each point in the phase diagram 
and thereby the transition lines between swinging (black), intermittent (grey shading) and tumbling regime (light grey) can be 
distinguished by the mean tumbling rate (8) defined by eqns. (|10I12|I which is proportional to the color shading. In panel a) the 
special line A = and the special point at A = and \ l — 1 which are discussed in section [3] are depicted explicitly. 



termittent" regime in Ref. [3] (see Fig. 2]). For low shear 
rates, i.e. for large values of the hydrodynamic flow 
is too weak to overcome the elastic barrier Eo due to the 
shape memory. Thus, the capsule tumbles for large x _1 - 
For smaller values of x _1 j i- e - f° r higher shear rates, the 
transition to the intermittent or swinging regime occurs. 

The dynamics in the reduced model can conveniently 
be characterized by investigating the normalized mean 
tumbling rate 



0) 



(&tB) 



(d T 9) + (d T <t>) 



(10) 



as an order parameter. Here the mean rates of inclination 
and phase angle are given by 

T 

(d T 6) = hm ^ J dTd T 6{ T ) and (11) 
o 

T 

(d T 4>) = r limi J dTd T 4>{ T ) , (12) 
o 

respectively. In a stable tumbling motion, the inclination 
angle grows without bounds while the phase angle oscil- 
lates, which implies (d T (f>) / (d T 9) = and a mean tum- 
bling rate (0) = 1 in the long time limit. Conversely, in 
a stable swinging motion, the phase angle grows without 
bounds while the inclination angle oscillates, which im- 
plies (d T 9) I (d T <f)) = in the long time limit thus (6) = 0. 
In the intermittent regime, the mean tumbling rate takes 
values between and 1. Typical phase diagrams showing 
grey scale plots of the mean tumbling rate (6) as obtained 
from solving equations JHJ and (J9j) numerically can be seen 
in Fig. [U Here, the axes consist of the dimensionless vis- 
cosity ratio A and the inverse capillary number x _1 j while 
the eccentricity a is constant. 

Even though it is not central for this paper, we note for 
completeness that the status of the intermittent regime is 



still under debatte. Kessler et al. [l8| and Sui et al. [20( 
solved the full dynamics of a 3d elastic capsule using a 
spectral method and a immersed boundary lattice Boltz- 
mann method, respectively. While the reduced model cap- 
tures the swinging and tumbling regime semi-quantitatively 
compared to fully numerical studies fl8l ] , the intermittent 
regime has been contested as an artifact of the reduced 
model. Neither study found any indications of intermit- 
tency, but rather a transition towards swinging. There was 
also no direct evidence of intermittent motion in experi- 
ments [1J. 



2.3 Quasi-spherical case 

Since the phase diagram is qualitatively similiar for all 
small values of a < 1 (see Fig. [U , it is sufficient to investi- 
gate the quasi-spherical case, for which analytical progress 
becomes possible. We set 0,2 = (1 — e)a\ where e <C 1 is 
a small parameter and assume the difference of 03 — a% 
to be also of order e. The three dimensionless parameters 
introduced above then depend on e to first order as (see 
Appendix [B] for definition of symbols) 



a rs e , 

3 + 277777° 



A 



X 



5Vrf 



-7E. 



(13) 
(14) 

(15) 



In the quasi-spherical case, the elastic energy can be cal- 
culated for any elastic model. In the regime of small de- 
formations the elastic energy scales quadratic with the 
eccentricity 



with 6q ~ 1, leading to 



(16) 



(17) 
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For given values of all physical parameters, leaving aside 
the shear rate 7 and the small parameter e, the pre- factors 
in the above expansions (fT3|) , (fT4|) , and (fP7)) are of the or- 
der of unity in the quasi-spherical limit e — > 0. In this case, 
we are restricted to a small left hand stripe in the phase 
diagram (see Fig. [4]) defined by A ~ e. Here, the disputed 
intermittent regime has no influence on the dynamics and 
can be ignored. Since the transition between tumbling and 
swinging takes place at x 1 ~ 1, we will later specify the 
shear rate to be of the order of the expansion parameter 
7 ~ e. For the following expansion we merely require the 
scaling of x -1 not to be smaller than e 1 . 

Before we expand the equations of motion (jHJ and ((9j 
in e, we introduce a suitable stretched, rotated, and trans- 
lated frame in the #</>-plane (see Fig. [2j with coordinates 



27 = 
A = 



-(*+•+!) 



7T 

4 ' 



(18) 
(19) 



where 27 is, up to a constant, the angle of a tracer particle 
with respect to the direction of the shear flow. In these 
coordinates, the mean tumbling rate (|10[) can be written 

as 



(d T A) + (8 T E) 

2 (d r 27) 



1 



(20) 



Finally, expansion of the equations of motion ((SJ) and §§§ 
up to first order in e (note that A ~ e, a ~ e) yields the 
quasi-spherical equations of motion in the new coordinates 

<9 T 27 = A, (21) 
d T A = 4sin27sin,4 + A + 2(x _1 - 1) cos (27- A). (22) 

With the initial condition Eq = 27(0) the solution of the 
first equation is 



27(r) - 27 + Ar , 



(23) 



i.e. a tracer particle moves with constant angular velocity 
with respect to the dimensionless time r. Assuming that 
the shear rate does not change sign, we can use the angle 27 
as a time quantity to arrive at an autonomous differential 
equation on the torus 

A— =4sin27sinZ\ + A + 2(x~ 1 - 1) cos (27 - A) , (24) 
a27 

where x _1 is now considered to be a function of 27 rather 
than t. Here, 27 and A are not restriced to the interval 
[0, 7r) but can take all real values, thereby accounting for 
the number of revolutions on the torus. 



3 Constant shear rate 

We first summarize the big picture derived below in de- 
tail for the shape dynamics expressed by 27 and A close 
to the special line A = (see Fig. |4j, corresponding to 



the quasi-spherical limit e — > at fixed inverse capillary 
number x ■ F° r a systematic expansion in A, we first in- 
vestigate the dynamics at the special line. We find a set 
of fixed points connected by curves, of which some are 
stable, some unstable, and some neutral (see Fig. [5]). For 
capillary numbers x 7^ 1) closed and separated curves in 
the phase space exist which consist solely of either stable 
or unstable fixed points. Upon perturbation to first order 
in A these lines turn to stable or unstable limit cycles, 
as the perturbation does not alter the absolute stability. 
These limit cycles correspond to swinging and tumbling 
for x 1 < 1 an d > 1) respectively. At the special 
point x _1 = 1, the lines of fixed points cross. Here, all 
vertical lines of fixed points consist solely of neutral fixed 
points, whereas all horizontal lines of fixed points con- 
sist of segments of either stable or unstable fixed points. 
Upon perturbation, the neutral fixed points can become 
either stable or unstable. Therefore the dynamics of the 
system close to the special point needs to be studied more 
carefully to first order in A. Analytic determination of the 
resulting limit cycles becomes possible by considering the 
trajectories close to the stable, unstable, and neutral fixed 
points separately and joining these with the method of 
matched asymptotic expansion. 



3.1 Zero order expansion on special line 

We start by investigating the special line A = 0, where the 
equations of motion 

<9 r 27 = 0, (25) 
d T A = 4 sin 27 sin Zi + 2(x _1 - 1) cos (27- A) (26) 

immediately lead to time-constant 27. For all values of x 
there are two connected lines of fixed points (see Fig. [5J 
which bifurcate into limit cycles upon a perturbation with 
A > 0. The position and stability character of these lines of 
fixed points depend on the inverse capillary number x _1 - 
The corresponding regimes on the special line A = are 
seperated by the special point (x _1 = 1, A = 0). 

At this special point, there are two straight lines of 
connected fixed points, namely 27 — kit with arbitrary A 
and A = kn with arbitrary 27 and integer k. These lines 
of fixed points build up a checkerboard pattern as can be 
seen in Figs. [5] c) and El We now discuss the stability of 
each fixed point. The eigenvalue of a linearization around 
a fixed point in the direction of the connected fixed points 
is always zero. The sign of the other eigenvalue determines 
the stability of the fixed point in the perpendicular direc- 
tion. Fixed points with a positive, negative or zero eigen- 
value are called unstable, stable or neutral, respectively. 
The vertical lines of fixed points are neutral as 27 is con- 
stant ([25]) . while stable and unstable segments alternate 
on the horizontal lines as can be seen in Fig. [5] c) , where 
the vector field corresponding to the equations of motion 
is shown. Thus, each square consists of two opposite neu- 
tral sides and two opposite lines of which one is stable and 
one is unstable. 
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For arbitrary points on the special line A = with 
X~ 1 7^ 1, there are two separated curves of fixed points 
(see Fig. [5] a, b) separated from the checkerboard pattern 

,i/2 

by a distance of order \x~ 1 ~ 1 as can be calculated 
using eqn. (|26p . Here, the fixed points on one single curve 
are either all stable or unstable. 



3.2 First order expansion away from the transition 

Since the system moves with constant velocity A along the 
A7-direction (see eqn. ([21])) . fixed points exist only for van- 
ishing viscosity contrast A = 0. For finite A > 0, the lines 
of fixed points turn into limit cycles. Excluding the region 
close to the special point (i.e. excluding x" 1 — 1 ~ A), the 
perturbation due to finite A > is too small to change 
the stability qualitatively. Thus, the stable character of 
the original stable line of fixed points as well as its topol- 
ogy remain unchanged. The stable limit cycle for small 
A > and x 1 < 1 leads to a decreasing A, while E is 
increasing (see Fig. [5] a). In the long time limit the mean 
rates have the same magnitude (d T A) = — (d T S), result- 
ing in a swinging motion with vanishing mean tumbling 
rate (9) = (sec eqn. ([20])). For > 1 the stable limit 
cycle leads to (d T A) = (d T S) > (see Fig. b), result- 
ing in a tumbling motion with mean tumbling rate (9) = 1 
(see eqn. ([20])). This is consistent with the phase diagrams 
shown in Fig. [U 



3.3 Matched asymptotic expansion close to the 
transition 

We now investigate the system for a finite but small vis- 
cosity contrast A ~ e close to the special point (x = 1, 
A = 0) , where we specialize to straight lines emerging from 
the special point. We define a slope parameter p ~ 1 of 
order unity by 



X 
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-A 



(27) 
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Fig. 6. The four regions i) to iv) with linear dimension of 
order vA in which asymptotic solutions are obtained in sec- 
tion (33] for < A < 1. i) Light grey region: thick stable line, 
ii) Dark grey region: junction region, crossing of stable and 
neutral line, iii) Light grey region: neutral line, iv) Dark grey 
region: crossing of neutral and stable line. 



(28) 



and solve the autonomous equation of motion ([24 
dA 4 



dS A 



• sin U sin A 



+p+{p-l)(cos(S-A) - 1) 

asymptotically. Since the curves of fixed points are sep- 
arated by the lines of fixed points of the checkerboard 
pattern by approximately a/Ix -1 ~ M ~ v^A, the stable 
limit cycle should be within stripes of width VA from the 
lines of fixed points of the special point. In each of the 
four regions i) to iv) shown in Fig. [SI we expand the equa- 
tions of motion up to lowest order in A ~ e and solve them 
analytically. By the method of asymptotic matching (36| 
the constants of integration can be deduced step by step. 
Here, we proceed as follows: We start with the general 
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i) ii) iii) 




-1 -0.5 0.5 -1 -0.5 0.5 -1 -0.5 0.5 



E/tt S/tt E/jr 

Fig. 7. i) to iv) Asymptotically matched solutions in the four regions close to the stable and neutral lines, v) All asymptotic 
solutions combined, vi) Numerical solution obtained by direct integration of the full equations of motion. Parameters are chosen 
from the tumbling regime (black curves) with a = 0.0017T, A = 0.3, p = 1 and from the swinging regime (grey curves) with 
p = — 1. The arrows denote the direction in time. 



solution of region i) and match the solution of region ii). 
We continue by matching region ii) with region iii) and 
region iii) with region iv). Finally, we match the solution 
of region iv) with the general solution of region i) to ob- 
tain a closed trajectory on the torus. Thus, we arrive at 
the unique analytical solution of the stable limit cycle. 

We now give a brief description of the solution which 
should be read while comparing Fig. [3 which shows typ- 
ical graphs of the obtained solutions. For a more detailed 
derivation, which includes the expanded equations of mo- 
tions and the asymptotically matched solutions, we refer 
the reader to Appendix ICl 

i) We start with the region close to the stable line with 
S ~ 1 and A ~ VA (see Fig. As is shown in Ap- 
pendix [Cl in region i) the limit cycle is to first or- 
der in A simply given by the original stable line (see 
eqn. (|C.3[) and Fig.[7]i). Thus, the system runs on the 
stable line A(U) — with — n < £ < irrespective 
of the slope parameter p. When the system starts in 
the vicinity of the limit cycle, it relaxes quickly to the 
stable line. 



ii) In the region S <~ vA and A ~ where stable 
and neutral line meet, the vector field corresponding 
to A = is small enough for the finite but small value 
of A to have a significant influence on the vector field 
and thus on the motion. Here, the exact value of the 
shear rate or slope parameter p is critical as can be 
seen by the matched solution (IC.5j) whose sign in the 
long-time limit depends only on the sign of the slope p 

Urn sign(^\(I7)) =signp. (29) 

A typical graph can be seen in Fig. [7]ii) for both cases 
p < and p > 0. For a negative slope p < the 
neutral line with A < is choosen which leads to a 
swinging motion. Conversely, for a positive slope p > 
the neutral line with A > is choosen, corresponding 
to a tumbling motion (see region hi)). Region ii) with 
S ~ A ~ y/X acts as a junction which only depends 
on the sign of the slope parameter p. 

iii) In the region U ~ VA and A ~ 1 close to the neutral 
line, the matched solution (|C.8|) describes the relax- 
ation towards the next stable line (see Fig. [7] iii), which 
has been chosen in region ii). 
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iv) In the region S ~ \f\ and A ± tt ~ \/A, where neutral 
and unstable lines meet again, the matched solution 
(|C.10P describes a relaxation toward the stable line for 
all values of p (see Fig. [7]iv). The system then starts 
over again in region i) close to the stable line. 

Panels v) and vi) of Fig. [7] show a comparison of the 
matched asymptotic solutions and the numerically com- 
putated stable limit cycle, with excellent agreement. 

Summarizing the dynamics, the system starts running 
along a horizontal stable line A = 0. At its end E ~ 0, 
it chooses one side depending on the sign of p and runs 
close to the vertical neutral line towards the neighbour- 
ing horizontal stable line. For negative p < 0, the angle 
A is decreasing along the neutral line. This case corre- 
sponds to a motion with oscillating inclination angle 9 
and monotonously decreasing phase angle <j) (see Fig. [7] 
vi), resulting in a swinging motion. For positive p > 0, 
the angle A is increasing along the neutral line. This case 
corresponds to a motion with oscillating phase angle <f> and 
monotonously decreasing inclination angle 9, resulting in 
a tumbling motion. 

In summary, these results imply for the phase diagram 
that the boundary between the tumbling and swinging 
regime is given by the line x" 1 = 1 — A/2, which corre- 
sponds to the critical value p = 0, in first order in e. 



4 Time-modulated shear rate 

We now investigate the dynamics of a capsule in a time- 
modulated shear flow and specialize to a periodically os- 
cillating positive shear rate ^(t) > with period T, fre- 
quency to = 2tt/T and mean value 



We first show numerical results for the mean tumbling 
rate (9) as a function of the driving frequency Co. These 
results were obtained by a direct integration of the equa- 
tions of motion (JS]) and (J9]). Fig.[8]shows the characteristic 
dependence on the frequency for a harmonically oscillating 
inverse shear rate. There is a smooth background, which is 
constant at low frequencies and vanishes at a high cut-off 
frequency. A large number of regularly ordered resonance 
peaks are superimposed. A qualitative discussion based 
on the results for time-constant shear rate can explain the 
general features of this plot. 



4.1 Qualitative explanation 

We start with some preliminary considerations, which will 
be confirmed analytically afterwards. As shown in sec- 
tion 13.31 for constant shear rates in the quasi-spherical 
limit, the equations of motion in regions i) and iii) and 
the qualitative relaxation towards the stable line in region 
iv) are independent of the shear rate. This behaviour re- 
mains unchanged for a time-dependent shear rate or slope 
parameter p{S). Thus, the relaxation towards the stable 
line and the motion on the stable line are unaffected by 
the shear rate. As illustrated in Fig. [9j the system there- 
fore runs on the torus with monotonously increasing angle 
E. It moves close to the stable line into the junction re- 
gion. Here it turns to one vertical side depending on the 
value of p and reaches a neighbouring horizontal stable 
line. During the motion the junction region is visited over 
and over again periodically in time r. We can label the 
junction with angles (Sk, Ak), where 



kit 



(34) 



7o= \ I "(/)'//. 



(30) 



Then the slope parameter p{S) defined by equations ([5|), 
(|2~T|) . and (|2"Tj) is periodic and can be written as 



p(S) = pa + pib(E) . 



(31) 



Here, the oscillatory function b(S) has a vanishing mean 
value 



7r/(I) 

J b(£)dE = 0, 



(32) 



is normalized to the maximum value max|6(i7)| = 1, and 

is periodic b(S) = b(S + tt/lo) with the dimensionless 
frequency 



470 



(33) 



Thus, a frequency of Co = 1 corresponds to a full rotation 
on the torus in Z'-direction. 



with integer index k counts the number of visits (see Fig. [9]) 
and Ak is an integer multiple of n counting the difference 
of the number of tumbling and the number of swinging 
motions. Starting with index k = 0, the system reaches 
the junction at consecutive angles Sq, E%, £2, ■ ■ ■■ 

The only difference to the time-constant case of sec- 
tion [3] happens in the junction region ii), where the sys- 
tem leaves the stable line to follow the neutral line. Here, 
the value of the instantenous shear rate determines for 
the overall behaviour. Since the shear rate is now time- 
dependent, the slope parameter p can take different signs 
each time the system is in the junction region and can even 
change signs several times within the junction region. 

We first want to estimate the time the systems spends 
in the junction region ii) and consider corresponding limit 
cases of the driving frequency Co. Measuring time with re- 
spect to the nondimensional time r, the speed of S is A. 
The time the system needs to return to the junction is 
of order 1/A. Since the junction region has linear dimen- 
sion of the order a/A, the time the system stays within the 
junction region is of the order -\/A/A = l/VA- Thus, the 
fraction of time the system is within the junction region is 
given by the order of V~\ and the corresponding frequency 
is of order 1 / VA- 
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a) b) 




0.2 0.4 0.6 0.8 1 2.5 5 7.5 10 12.5 15 17.5 



Fig. 8. Typical numerical value of the mean tumbling rate (0) as a function of <I> in finite simulation runs of time-dependent shear 
flow, a) Low frequencies < Q < 1. b) Whole spectrum. - There are several resonance peaks of different height whose width 
depends upon both the length of the simulation run and the amplitude of the oscillating shear rate. The smooth background 
starts at a constant value for small frequencies and finally vanishes at a cut-off frequency. The thick grey line is an analytical 
result (|44[) describing the background with excellent agreement. Parameters (see eqns. (|31|l . (138[) and sec. I4.3|l : a — 0.0017T, 
A = 0.01, po = -1, pi = 2, cp = 0, E f = 500tt. 



For high frequencies uj 3> 1/ \/\, the oscillation is too 
fast for the system to respond. Therefore the shear rate be- 
haves effectively as a time-constant shear rate with mean 
slope po. For a negative mean slope po < 0, there is a 
pure swinging motion with vanishing mean tumbling rate 
(6) = (see Fig. [8]b). Conversely, for a positive mean 
slope pq > there is a pure tumbling motion with mean 
tumbling rate (9) = 1. 

In the limit of low frequencies Q <C 1/ \f\, the shear 
rate in the junction region can be regarded constant. In 
other words, the junction region is effectively just a point 
located at (Sk, Ak). Each time the system is in the junc- 
tion region labeled by Sf. , the sign of the slope parameter 
p(Sk) determines whether the system performs a single 
tumbling or a single swinging motion. The sign of p(SJk) 
depends on the initial phase ipo, the frequency to, and 
the index k. In order to calculate the mean tumbling rate 
(6), the number of positive and negative values of p(So), 
p{S\), . . . have to be counted. Therefore, the system can 
be mapped on a discrete model in the low-frequency limit 
as shown in section |4~3")) . This discrete model reproduces 
both the constant background of the mean tumbling rate 
and the superimposed resonance peaks which can both be 
seen in Fig. [5] 

For intermediate frequencies Cj ~ 1/V% it seems rea- 
sonable that some time-averaged slope p in the junction 
region determines the motion of the capsule. This expec- 
tation is quantified in the next section. 



4.2 Analytic solution in the junction region 

The qualitative arguments of the previous section are sub- 
stantiated by a full analytical investigation for a general 
time-dependent shear rate with parameters A and x _1 
close to the special point. Since we are in the quasi-spherical 



<T 



I I 
I I 

<(2 lT A r j/7r- 



V I- 



S 4 ,A 4 )/7r 

i 



1 2 3 4 5 6 

E/tt 

Fig. 9. Sequence of junctions labeled by (-£"*,, At) f° r a given 
trajectory. 



limit, we investigate the four regions analog to the time- 
constant case (sec. l3.3j) . The equations of motion to lowest 
order in A remain unchanged except for the fact that the 
slope parameter is now time-dependent p = p{£). Since 
the velocity in I7-direction in dimensionless units is A, the 
system visits the junction labeled by (£k, ^fe) at time 
t = irk/X. 

As in the time-constant case (see Appendix [C|l . there 
is no dependence of the solutions on the slope parameter 
p in regions i) and iii), and no qualitative dependence on 
p in region iv), where the trajectory merely relaxes to the 
next stable line. Thus, after leaving the junction region, 
the system moves fast towards one of the two neighbour- 
ing stable lines, before moving slowly along the stable line 
and returning to the junction region ii). Again, the junc- 
tion region determines which stable line is chosen next, 
i.e. whether the capsule tumbles or swings. The corre- 
sponding first order equation of motion (|C.4[) in region ii) 
close to Sk can be integrated for a general time-dependent 



10 



Steffen Kessler et al.: Elastic capsules in shear flow 



p{S) as is shown in Appendix [D] In the long-time limit 
r ^> l/VX the matched solution (|D.2|) in the junction 
region becomes asymptotically 



A(E) - A k 



exp 



s h ) 2 



Pk 



(35) 



where the average slope pk corresponding to the k-th junc- 
tion (Skj Ak) is defined by 



Pk 



i 00 

dEp(Z)exv f-l{S-S k ) 2 \ 



(36) 



This integral is a convolution of the time-dependent slope 
parameter p{S) with a Gaussian shaped kernel of width 
vA centered at S = Sk- Thus, for low frequencies Co <C 
1/VA the kernel is effectively proportional to Dirac's 8- 
function, while for high frequencies Q 3> 1/vA the kernel 
smoothes out the fast oscillations of p{S). These two limit 
cases will be discussed in more detail in the next section. 

The sign of the average slope pk determines whether 
the trajectory of the system follows the neutral vertical 
line along the positive (for p~k > 0) or negative (for pk < 0) 
direction. Since the solutions in regions iii) and iv) only 
describe the relaxation to the next stable line, the asymp- 
totic matching procedure then proceeds exactly as in the 
time-constant case. For any given time-dependent shear 
rate ^(t) or equivalently p{£), the sequence of average 
slopes p~k at S — S k can be calculated. The mean tum- 
bling rate is then given by 



N-l 



(0) 



lim 

N-*oo N 



(37) 



fc=0 



with the Heaviside step function 0. We now evaluate this 
expression for a specific choice of p{S). 



4.3 Harmonically oscillating shear rate 

Since p{£) is periodic with period ft /Co, it can be de- 
composed into a Fourier series consisting of an oscillation 
with the fundamental frequency 2Co and the correspond- 
ing higher harmonics. We constrain p{£) to a pure har- 
monic oscillation in the following section for simplicity. 
The results are easily generalised to the Fourier series of 
an arbitrary periodic p{£) (see Appendix |E|) . 
For a purely harmonically modulated slope 



p{E) = p Q + pi cos (2(u)S + (p )) 



(38) 



with mean value po, amplitude p\ > 0, frequency Co and 
initial phase ipo, the sequence of mean slopes pk can be 
evaluated analytically 

Pk = Po + cxp J Pi cos(2ip k ) , (39) 



where we have defined the sequence of equidistant phases 
fk = fa + kiruj . (40) 
We now discuss the mean tumbling rate 

m 

N- 



(0) = lim (Po + ^P (~\Cj 2 /2) Pi cos(2^ fc )) (41) 



k=0 



as a function of po, p\, Cj, and (po . 

For small modulation amplitudes exp | 



A£ 2 /2) Pl 



< 

I Pol, the sign of po equals the sign of pk for all k. This 
condition holds for all frequencies if pi < \ P q\. For p\ > 
\ P o\, we can define a threshold frequency 



A" 1 / 2 4 2 In 



Pi_ 

\Po\ 



(42) 



beyond which {Co > Co c ) the phase behaviour is given by 
the sign of po alone. The system tumbles ((9) = 1) for 
po > and swings ((9) = 0) for p < 0. 

For large modulation amplitudes exp (— A(D 2 /2) p\ > 
\ P o\, the sign of pk depends on the value of the phase <pk 
modulo 7T. For phases in the region \ifk\ < Aip/2 with 



A(p 



arccos 



( p° 1 






:")) 





(43) 



the average slopes pk are positive and the system performs 
tumbling motions and swinging motions otherwise. This 
behaviour is visualised in Fig. 1101 where the phases ifk 
modulo 7r are interpreted as points exp(2iipk) on the circle 
S 1 (modulo 7r). 

The circle S 1 consists of an arc with angle A<p corre- 
sponding to tumbling and a complementary arc with angle 
7r — A(p corresponding to swinging. Each phase tfk at the 
junction labeled by {Sk, Ak) is either an element of the 
tumbling or an element of the swinging arc. By counting 
the fraction of phases within each arc the mean tumbling 
rate can be calculated explicitly: 

— For an irrational frequency Co the values of the phases 
= Po+kirCo lie densely on the circle S . The fraction 
of number of swinging to tumbling motions in the long 
time limit is then given by the ratio of the length A(p 
and 7r — Aip of the two intervals, leading to a mean 
tumbling rate of 



(0) 



Aip 

7T 



arccos (— po exp(Aa) 2 /2)/pi) 



(44) 



In the low frequency limit this becomes a constant 
(6) « arccos(-po/pi)/7r- 

For a rational frequency Co = p/q with integer and co- 
prime numbers p and q the phases ifk — + knCo — 
tpo + kpn/q lie on q equidistant phases ipo +jTt/q (with 
j = 0, . . . , q — 1) on the circle S 1 . In the long time 
limit, all of these q angles are visited equal amounts 
of times. The ratio of number of swinging to tumbling 
motions is given by the ratio of number phases ipk in 
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the two intervals of length Atp and tt — Aip. For high 
values of q, approximately qAip phases lie within the 
interval Aip. Then the mean tumbling rate (9) is ap- 
proximately given by the expression (|44p valid for irra- 
tional frequencies. For low denominators q, the number 
of phases within the two intervals additonally depends 
upon the initial phase ipo . An integer frequency Q = p 
for instance gives a tumbling rate of either (8) = 1 or 
(9) = depending only on the initial condition. For 
a general rational frequency, counting the number of 
phases \ipk\ < Atp in the tumbling sector gives 



0) = 



1 



Aip 

TT \~Y~ 

' (45) 

with Aip given by eqn. (|43|) . Here, |_ J and [ ] are the 
floor and ceiling functions, which denote the closest in- 
teger smaller or larger than the argument, respectively. 

Plotting the mean tumbling rate {6} over small frequen- 
cies Q in Fig. [8j we can identify a smooth irrational back- 
ground (|44p superimposed by rational peaks (|45|). Their 
amplitudes depends on the denominator q and the initial 
phase ipo. For integer resonance frequencies, the peaks go 
either to or 1. 

Mathematically speaking, the peaks at rational fre- 
quencies are infinitesimally narrow. Since experiments and 
numerical methods run only for a finite time, the peaks 
are broadened to a finite width, which can be estimated 
in the following for the most dominant peaks with nearly 
integer frequency 

Q = p + SQ. (46) 

Let us assume that ipo = 0, so that the system tumbles 
with (8) = 1 (for Q = p). For exactly integer modulation 
frequency SQ = 0, the phase stays constant ipk = for 
all fc. For a non-integer frequency SQ ^ the fc-th phase 
modulo 7r is pu = kirSQ = SQUk- If this phase does not 
change by more than by the width of the tumbling arc Atp 
during the simulation time Sf, the system stays in the 
tumbling regime. Therefore we can estimate the width of 
the integer resonance peaks as 



AQ 



Aip 



1 



( p° 1 












b) 




(47) 



Fig. 11. a) Mean tumbling rate (9) at the rational peak cor- 
responding to Q = 1 for different slope amplitudes p\ = 2., 
1.4, 1.1, 1.04. The higher the amplitude pi, the smaller the 
width Aui. The constant parameters are a = 0.0017T, A = 0.01, 
po = — 1, a finite running time Sf — 40tt, and the initial phase 
ipo within the tumbling arc. b) Numerical (dots, from a) with 
{&) > 0.999) and analytical (solid line, eqn. (|47p ) results of the 
peak widths Au) for different amplitudes p\ (remaining param- 
eters as in a). 



This result is confirmed in Fig. II H where numerically eval- 
uated peak widths are compared to expression (J47J) for dif- 
ferent modulation amplitudes p\ . For the numerical data, 
the peak width was defined to be the frequency interval in 
which the mean tumbling rate during the simulation time 
(0) was larger than 0.999. The agreement is excellent. 
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4.4 Dynamic phase diagram 

We summarise our findings in a dynamic phase diagram. 
While the specific shape of the phase diagram depends 
on the particular functional time-dependence of the shear 
flow, the general features hold for any oscillating time- 
dependent flow with mean inverse shear rate Xq 1 an( i 
oscillation amplitude Xx i compare cqn. (|31j) . From our 
study of time-constant shear flow we know the location 
of the phase boundary x^ 1 = 1 — A/ 2. In Fig. [121 a grey 
scale plot of the mean tumbling rate (9) is shown as a 
function of the oscillation amplitude Xi relative to the 
distance of the mean shear rate to the phase boundary 
Xi 1 /\Xq 1 ~~ Xq U versus the oscillation frequency ui — 
Lv/4rfQ measured im units of the mean shear rate. The 
colour level indicates the mean tumbling rate as defined 
in eqn. (jTUJ) for finite simulation times, black colour indi- 
cating a swinging motion (9) — 0. In Fig.[l2j a mean shear 
rate in the swinging regime (for constant flow) was chosen, 
X^ 1 < X7 1 - One can see that for oscillation amplitudes be- 



low the distance to the phase boundary, X\ 1 < Xc 1 ' 



Xo 1 



the capsule never tumbles. In order to induce tumbling 
motion, the instantaneous shear rate has to cross the phase 
boundary. For higher modulation frequencies, the oscilla- 
tion amplitude threshold for tumbling is even higher and 
given by 



Xil 



Xo 1 ! 



exp 



Xu 2 
32 7 2 



(48) 



Above the oscillation amplitude threshold the mean tum- 
bling rate grows continuously on the irrational background 
with increasing amplitude. At the resonance frequencies 
the mean tumbling rate reaches values given by eqn. (|45l) . 



5 Conclusions 

We have investigated the motion of microcapsules in time 
dependent shear flow in a reduced model. The equations 
of motions were studied analytically in the quasi-spherical 
limit for constant viscosity contrast. We have identified 
the stable and unstable fixed points at lowest order in the 
deformation, which lead to swinging and tumbling limit 
cycles at first order depending on the shear rate. Close to 
the dynamic phase boundary, the expansion of the equa- 
tions of motion was carried out to first order in the defor- 
mation. Their analytic solution was facilitated by solving 
the trajectories in different regions in phase space sepa- 
rately: One region close to the stable fixed points with 
comparatively slow dynamics, one region close to the un- 
stable fixed points with a comparatively fast motion, and 
a junction close to the neutral fixed point, which also acts 
on a fast time scale. The direction of the dynamics in this 
latter region is determined by the value of the shear rate. 
The analytic inner solutions for the trajectories were then 
joined together using the method of matched asymptotic 
expansion. 



X 



I u 
X 



X 




6 8 10 
w/(4 7 ) 
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Fig. 12. Dynamic phase diagram visualized as a grey scale 
plot of the mean tumbling rate {&) as a function of oscillation 
amplitude relative to the distance of the mean shear rate to the 
phase boundary xT 1 / (XtT 1 — Xo" 1 1 > an d 01 the normalised mod- 
ulation frequency w/(4<yo). Black colour indicates pure swing- 
ing, (6) = 0. Tumbling only occurs above a threshold in the 
oscillation amplitude given by eqn. (|48p . Data shown in Fig. [5] 
correspond to a cut in parameter space as indicated by the 
grey horizontal line. 



As a central result of this study, we now fully under- 
stand the dynamic phase behaviour of quasi-spherical cap- 
sules in time-independent shear flow and have determined 
the phase boundary between swinging and tumbling as a 
function of shear rate and viscosity as = 1 — A/2. In 
physical parameters the phase boundary reads 



7 



Vf]°£ 



10 



(49) 



We then generalised our result to time-dependent shear 
rates. The equations of motion were integrated analyti- 
cally in the different dynamic regions. Again, we found 
that the motion for one period is determined by a weighted 
time-average of the shear rate during the time when the 
system is close to the junction. We have thus mapped the 
continuous capsule dynamics to a discrete model valid for 
all time-dependent flows. Our general results were then 
applied to harmonically modulated shear rates around 
a finite mean value, where the dynamic phase diagram 
was constructed explicitly. As the dynamics is determined 
by the (time-averaged) value of the shear rate at specific 
times only, the system shows a pronounced resonance be- 
haviour: By choosing suitable modulation frequencies, it 
is possible to induce tumbling motion for capsules, which 
would otherwise swing at constant mean shear rate. The 
width of the resonance peaks for finite simulation time can 
also be expressed analytically. For not resonant frequen- 
cies, the mean tumbling rate is determined by an analytic 
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expression, which vanishes beyond an upper modulation 
frequency threshold. The agreement of numerical solutions 
of the original equations of motion with the theoretical re- 
sults is excellent. 

In summary, we have reached a fairly complete ana- 
lytical understanding of the motion of a quasi-spherical 
capsule in time-dependent shear flow. In the course of 
the study, we have also obtained analytical expressions for 
capsules in constant shear flow whose equations of motion 
previously have been studied merely numerically. 
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A Keller-Skalak quantities 

The geometrical quantities /, occuring in equations (|2j) 
and ^ depend upon the semi-axes oi, 02 1 03 and are 
defined by the following relations: 



d2 a 3 
r-i = — , r 3 = — , an = 

°i °i (010203) 



A = v /(a? + s)(ai + s)(a§ + s), 

00 

ds 

93 



(a\ + s)(a| + s)A 



(A.2) 
(A.3) 

(A.4) 



2r 2 2 

fx = (r 2 - — V . h = ^U--) (A.r.) 



r 2 



z-2 



Z2 



(A.6) 



In the axisymmetric cases a 2 — 03 or a\ — 03 and in the 
quasi-spherical case |1 — 02,3/01! £ <C 1, the integral g' 3 
can be computed explicitly. 



B Notation for asymptotic limits 

We use the following conventions, where f(x) and g(x) are 
real functions of the real variable x, c is a real constant 
and y and z are real numbers: 



f(x)~g(x) O lim4^ = <. 



x^Q g(x) 

m 



1 



(B.l) 



(B.2) 



f{x) ~ g(x) & lim 

x^o g(x) 

y ~ z y and z numerically equal . (B.3) 



special point A = and = 1- For each of the four 
regions of Fig. [5J the quasi-spherical equation of motion 
(|29|) is solved asymptotically, compare Fig. [7j 



i) S ~ 1, A ~ VA: Expansion gives the equation of 
motion 



dA 4 A 

ds = x As[nI] 



(C.l) 



which is independent of the shear rate, i.e. independent 
of the slope p. Its solution is given by 



A{E) = Ai exp (^-^(cos S - cos Si) j , (C.2) 

where (Si, A{) is an arbitrary point on the trajectory. 
This solution will be used to match with the solution 
of region ii). As can be seen in iv) by closing the trajec- 
tory on the torus, the limit cycle in region i) is simply 
given by the original stable line 



A(S) = 0. 



(C.3) 



Thus, to leading order in A ~ e, the system runs on 
the stable line A(S) = with -n < S < 0. Even 
when the system starts off the limit cycle, the stable 
character of A — and — ir < S < leads to a fast 
relaxation towards the stable line, while the angle S 
changes slowly due to A < 1. Therefore, the match- 
ing with region ii) will not depend on the initial point 
(Si, Ai) and the general solution (|C2[) is independent 
of the slope parameter p. 
ii) S ~ VA, A \f\: Here, stable and neutral line meet, 
and the expansion gives 



4 . ^ 

ds = x AS + p - 



(C.4) 



In this region, the vector field corresponding to A = 
is small enough for the finite but small value of A to 
have a significant influence on the the vector field and 
thus on the motion. Here, the exact value of the the 
shear rate or the slope parameter p is critical as can 
be seen by the solution 

/ ttA 

A(E) = \l—p 



(C.5) 




x exp 



which has been matched with the general solution (|C.2[) 
of region i) and which is independent of the initial 
point. Here, we use the error function 



C Asymptotic matching for time-constant 
shear rates 

Here, we describe our analytical procedure for obtaining 
the asymptotic trajectory of the limit cycle in the quasi- 
spherical case A ~ e -C 1 when the system is close to the 



X 

evtx=—= I dsexp(—s 2 ). 
V" J 



(C.6) 



Thus, the sign of p determines the sign of A(E). For 
negative p < the neutral line with A < is choosen 
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which leads to a swinging motion. Conversely, for a 
positive p > the neutral line with A > is choosen, 
corresponding to a tumbling motion (see region iii)). 
Region ii) with E ~ A ~ %/A acts as a junction which 
only depends on the sign of the slope parameter p. 
iii) In the region E ~ vA and A ~ 1 close to the neutral 
line, the expansion gives a slope-independent equation 
of motion 



^ — oo has to be taken. The matched solution in the 
junction region then is given by 



A(E) = exp 



2 9 
A 



(D.2) 



dE'p{E') exp ( — — i7 



d/i 4 

= -A7smZ\, 

dA7 A 



whose matched solution 



A(E) = 2arctan ^tan(^/ — p) 
2 



(C.7) 



(C.8) 



For a harmonically changing slope parameter (|38|) with 
mean value po, amplitude pi, frequency tu and initial phase 
ipo the matched solution (|D.2j) can be integrated 



A(E) 



x exp 



depends on the slope parameter p. This solution de- 
scribes the relaxation towards the next stable line which 
has been chosen in region ii). 
iv) In the region E ~ y/\ and A ± tt ~ VA, where neutral 
and unstable lines meet again, the expansion gives 




(D.3) 



•2-p 



(C.9) 



cos 2ipo Re |^erfc 

sin 2<po Im (eric 

where erfc z = 1 — erf z is the complementary error func- 
tion, Re and Im denote real and imaginary parts, respec- 
tively. 



with matched solution 



A(E) = T7i" + 2exp --E 



±7r — cot(y— p) 




(CIO) E Fourier series 



(2-p)ern(^/-i:) 



where we used the imaginary error function 



erfix 



( T i [ix) 



(C.ll) 



The upper and lower signs depend on which neutral 
line was chosen in region ii). Although both the equa- 
tion of motion and the solution depend on p, the sys- 
tem relaxes toward the stable line for all values of p. 
It then starts over again in region i) close to the stable 
line. By matching with the general solution of region i) 
the trajectory closes and the solution (|C.3[) is obtained. 



The oscillating part b(E) of any general period slope p(E) 
(see eqn. pip) can be uniquely decomposed into a Fourier 
series 

oo 

b(E) = ^ f b j coa(2(ju)E + v j )). (E.l) 
i=i 

In this case, the sequence of mean slopes pk (see eqn. (|3l>)) ') 
is explicitly given by the Fourier series 

Pk = Po + Pi X! h J ex P ( o — ) cos ( 2 (^ fc7r + u j))- 

(E.2) 

For b\ — 1 and bj — for all j > 1, we recover eqn. (j3"5|) . 
The amplitude bj is damped by the factor exp (— Aj 2 cD 2 /2). 



D Asymptotic matching for time-dependent 
shear rates 

The first order equation of motion (|C.4p in region ii) 

^ = iAE + p(E) (D.l) 

can be integrated for a general time-dependent p(E). In 
order to match asymptotically with the solution A(E) = 
(|C.3p of region i), the limit r — > — oo or equivalently 
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